Dark Radiation and interacting scenarios 
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An extra dark radiation component can be present in tiie universe in the form of sterile neutrinos, 
axions or otiier very light degrees of freedom which may interact with the dark matter sector. We 
derive here the cosmological constraints on the dark radiation abundance, on its effective velocity 
and on its viscosity parameter from current data in dark radiation-dark matter coupled models. 
The cosmological bounds on the number of extra dark radiation species do not change significantly 
when considering interacting schemes. We also find that the constraints on the dark radiation 
effective velocity are degraded by an order of magnitude while the errors on the viscosity parameter 
are a factor of two larger when considering interacting scenarios. If future Cosmic Microwave 
Background data are analysed assuming a non interacting model but the dark radiation and the 
dark matter sectors interact in nature, the reconstructed values for the effective velocity and for the 
viscosity parameter will be shifted from their standard 1/3 expectation, namely = 0.34to'oo3 
and cli^ = 0.29lJii:li!J5i at 95% CL for the future COrE mission data. 

PACS numbers: 98.80.-k 95.85. Sz, 98.70.Vc, 98.80.Cq 

I. INTRODUCTION 

From observations of the Cosmic Microwave Back- 
ground (CMB) and large scale structure (LSS) we can 
probe the fundamental properties of the constituents of 
the cosmic dark radiation background. The energy den- 
sity of the total radiation component reads 
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where is the current energy density of the CMB and 
A^eff is a free parameter, defined as the effective number 
of relativistic degrees of freedom in the cosmic dark radi- 
ation background. In the standard scenario, the expected 
value is AeS = 3.046, corresponding to the three active 
neutrino contribution and considering effects related to 
non-instantaneous neutrino decoupling and QED finite 
temperature corrections to the plasma. The most recent 
CMB data analyses gives A^eff = 3.89±0.67 (68% CL) [J, 
see also Refs. |2h18|. The simplest scenario to explain the 
extra dark radiation AA^off = A'off — 3.046 arising from 
cosmological data analyses assumes the existence of ex- 
tra sterile neutrino species. However, there are other 
possibilities which are as well closely related to minimal 
extensions to the standard model of elementary particles, 
as axions, extra dimensions or asymmetric dark matter 
models. 

Dark radiation, apart from being parametrized by its 
effective number of relativistic degrees of freedom, A'cff, 
is also characterized by its clustering properties, i.e, its 
rest-frame speed of sound, c^jj, and its viscosity param- 
eter, c^jg, which controls the relationship between veloc- 
ity/metric shear and anisotropic stress in the dark ra- 
diation background A value of cjjg different from 
zero, as expected in the standard scenario, sustains the 



existence of dark radiation anisotropies [20|. The stan- 
dard value of c^jg = 1/3 implies that the anisotropies in 
the dark radiation background are present and they are 
identical to the neutrino viscosity. On the other hand, 
the case cjjg = cuts the Boltzmann hierarchy of the 
dark radiation perturbations at the quadrupole, repre- 
senting a perfect fluid with density and velocity (pres- 
sure) perturbations exclusively. A value of c^g different 
from the canonical c^g = 1/3 leads to a non-adiabatic 
dark radiation pressure perturbation, i.e. {Sp—Sp/3)/p = 
{c^g — l/3)(52^''*, where S^'^'^* is the density perturbation 
in the rest frame, where the dark radiation velocity per- 
turbation is zero. 

Interacting dark radiation arises naturally in the so- 
called asymmetric dark matter models (see e.g. [2l[ and 
references therein), in which the dark matter produc- 
tion mechanism is similar and related to the one in the 
baryonic sector. In these models, there exists a particle- 
antiparticle asymmetry at high temperatures in the dark 
matter sector. The thermally symmetric dark matter 
component will annihilate and decay into dark radiation 
degrees of freedom. Since the dark radiation and the 
dark matter fluids are interacting, there was an epoch 
in the early Universe in which these two dark fluids 
were strongly coupled. This results in a tightly coupled 
fluid with a pressure producing oscillations in the mat- 
ter power spectrum analogous to the acoustic oscillations 
in the baryon-photon fluid before the recombination era. 
Due to the presence of a dark radiation-dark matter in- 
teraction, the clustering properties of the dark radiation 
component can be modifled [lll |. In other words, if dark 
radiation is made of interacting particles, the values of 
the clustering parameters c^g and c^jg may differ from 
the canonical c^g = c^j^ = 1/3. 

The cosmological implications of interacting dark ra- 
diation with canonical clustering properties have been 
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carefully explored in Refs. [23 - l23 |. see also the recent 
work of Ref. |18|]. Here we generalize the analysis and 
leave the three dark radiation parameters AA'cff, c^g and 
c^ig to vary freely within a ACDM scenario with a dark 
radiation-dark matter interaction. We will see that the 
bounds from current cosmological data on the dark ra- 
diation properties derived in non interacting schemes in 
Refs. 11, 13, 17] will be, in general, relaxed, when an in- 
teraction between the dark radiation and the dark matter 
fluids is switched on. While the bounds on the number 
of extra dark radiation species will not be largely modi- 
fied in coupled schemes, the errors on the dark radiation 
effective velocity and viscosity parameters will be dras- 
tically increased in interacting scenarios. We also show 
here how future CMB measurements, as those from the 
Planck [25,] and COrE [26] experiments, can lead to large 
biases on the dark radiation clustering parameters if the 
dark radiation and dark matter fluids interact in nature 
but the data is analyzed in the absence of such a coupling. 

The paper is organized as follows. Section HIl presents 
the parametrization used for dark radiation, describing 
the dark radiation-dark matter interactions explored here 
and their impact on the cosmological observables used in 
the analysis, as the CMB temperature anisotropies and 
the matter power spectrum. In Sec. HIT] we describe the 
data sets used in the Monte Carlo analyses presented in 
Sec. lIVi which summarizes the constraints on interacting 
dark radiation properties from current cosmological data. 
Future CMB dark radiation measurements are presented 
in Sec. El We draw our conclusions in Sec. |Vll 

II. DARK RADIATION-DARK MATTER 
INTERACTION MODEL 

The evolution of the dark radiation linear perturba- 
tions reads [l^ 

Sdr - ^(1 - Scl^) (Sdr + 4^^) + -^e^r + ^ = ; 

(2) 
(3) 

ndr + - Icl,, (^^+a^^0; (4) 

21 + 1 ■ 

—^Fdr^S - lFdrJ-1 + + l)Fdr,l + l =0 / > 3 , (5) 

where the dots refer to derivatives with respect to confor- 
mal time, a is the scale factor, k is the wavenumber, c^jj 
is the effective sound speed, cjjg is the viscosity param- 
eter, Sdr and 9dr are the dark radiation energy density 
perturbation and velocity divergence, respectively, and 
Fdr^i are the higher order moments of the dark radia- 
tion distribution function. In the set of equations above, 
TTdr is the anisotropic stress perturbation, and a is the 



metric shear, defined a.s a = {h + 6r])/{2k), with h and r] 
the scalar metric perturbations in the synchronous gauge. 
The anisotropic stress would affect the density perturba- 
tions, as in the case of a real fluid, in which the stress 
represents the viscosity, damping the density perturba- 
tions. The relationship between the metric shear and the 
anisotropic stress can be parametrized through a "viscos- 
ity parameter" , cjj^ (19(1 ■ 

TT = -3^7r + 4^^ (e + a) , (6) 

where 6 is the divergence of the fluid velocity. Although 
the perturbed Einstein and energy-momentum conserva- 
tion equations are enough to describe the evolution of 
the cosmological perturbations of non-relativistic parti- 
cles, it is convenient to introduce the full distribution 
function in phase space to follow the perturbation evo- 
lution of relativistic particles, that is, to consider their 
Boltzmann equation. In order to determine the evolu- 
tion equation of dark radiation, the Boltzmann equation 
is transformed into an infinite hierarchy of moment equa- 
tions, that must be truncated at some maximum mul- 
tipole order Imax- Then, the higher order moments of 
the distribution function are truncated with appropriate 
boundary conditions, following Ref. p7| . 

In the presence of a dark radiation-dark matter inter- 
action, the Euler equations for these two dark fluids read 

Sdm = ^ '—Odm + -J- 0.ndm<^dm-dr{ddr — &dm) , (7) 

a oPdm 

Sdr = ^fc^ {Sdr — "^T^dr) + andm^^dm-drif^dm — (^drX^) 

where the momentum transfer to the dark radiation com- 
ponent is given by andmC^dm-driddm - ddr)- Indeed, the 
former quantity is the differential opacity, which gives the 
scattering rate of dark radiation by dark matter [l^, ■ 
The complete Euler equation for dark radiation, includ- 
ing the interaction term with the dark matter fluid, reads 

Q q;22/^1c , a9dr\ CI 1 2 

Odr = 3k Cpff -ddr H pr &dr - -^1^ TTrf,. 

\4 a k/ J a 2 

+ 0'ndm<^dm.-dr{(^dm — Sdr) ■ (9) 

Following Refs. [13, [2^ we parameterize the coupling be- 
tween massless neutrinos and dark matter through a cross 
section given by 

{crdm-dr\v\) ^ Qorridm , (10) 
if it is constant, or 

{crdm-dr\v\) ^ ^mdm , (H) 

if it is proportional to T^, where the parameters Qo and 
Q2 are constants in cm^ MeV~^ units. It has been shown 
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in Ref. [2J] that the cosmological imphcations of both 
constant and T-dependent interacting cross sections are 
very similar. Therefore, in the following, we focus on the 
constant cross section case, parameterized via Qo- 

Figure [U upper panel, shows the CMB temperature 
anisotropies for Qq = 10^"^^ cm^MeV^^ and one dark 
radiation interacting species, i.e. AA'ofr = 1, as well 
as for the non interacting case for the best fit param- 
eter values from WMAP seven year data analysis d, |33| 
together with WMAP and South Pole Telescope (SPT) 
data We illustrate the behavior of the temperature 
anisotropies for different assumptions of the dark radi- 
ation clustering parameters. Notice that the presence 
of a dark radiation-dark matter interaction enhances the 
height of the CMB peaks due to both the presence of 
an extra radiation component (AA^off) and the fact that 
dark matter is no longer pressureless (due to a non zero 
Qo). Therefore AA'cff and Qq will be negatively corre- 
lated. The location of the peaks also changes, mostly 
due to the presence of extra radiation AAeff- The peaks 
will be shifted to higher multipoles £ due to changes in 
the acoustic scale, given by 

9a = , (12) 

where rg{zrec) and rs{zrec) are the comoving angular di- 
ameter distance to the last scattering surface and the 
sound horizon at the recombination epoch Zrec, respec- 
tively. Although re{zrec) almost remains the same for 
different values of AA'off, rs{zrec) becomes smaller when 
AA'off is increased. Therefore, the positions of the acous- 
tic peaks are shifted to higher multipoles (smaller an- 
gular scales) if the value of AA^eS is increased. Notice, 
however, that this effect can be compensated by chang- 
ing the cold dark matter density, in such a way that Zrec 
remains fixed, see Ref. ^9|]. Changing c^jg modifies the 
ability of the dark radiation to free-stream out of the po- 
tential wells ^28l-l30| . Notice from Fig. [T] (upper panel) , 
that lowering c^-^ to the value c^j^ = 0, the TT power 
spectrum is enhanced with respect to the standard case 
without the dark radiation and the dark matter species 
interacting. This situation can be explained, roughly, as 
the dark radiation component becoming a perfect fluid. 
That is, we are dealing with a single fluid characterized 
by an effective viscosity. Disregarding the fluid nature 
and the physical origin of the viscosity, the general con- 
sideration holds: for a given perturbation induced in the 
fluid, the amplitude of the oscillations that the viscosity 
produces (see, e.g. |Tl|) increases as the viscosity is re- 
duced. Therefore, lowering c^jg diminishes the amount 
of damping induced by the dark radiation viscosity, and, 
consequently, in this case, the amplitude of the CMB os- 
cillations will increase, increasing in turn the amplitude 
of the angular power spectrum. Therefore, we expect the 
interaction strength size Qo and the c^jg parameter to be 
positively correlated. 

On the other hand, a change of c^jj implies a decrease 
of pressure perturbations for the dark radiation compo- 



nent in its rest frame. As shown in Fig. [T] (upper panel), a 
decrease in c^fj from its canonical 1/3 to the value clg = 
leads to a damping of the CMB peaks, since dark radia- 
tion is behaving as a pressureless fluid from the pertur- 
bation perspective. In the case of c^g, we expect this 
parameter to be negatively correlated with Qo- 

Figure [T] (lower panel) depicts the matter power spec- 
trum in the presence of a dark radiation-dark matter in- 
teraction for different values of the dark radiation cluster- 
ing parameters (including the standard case with Cgfj = 
c^jg — 1/3) for one dark radiation interacting species, i.e. 
AA'off = 1. We illustrate as well the case of a pure ACDM 
universe. Notice that, since the dark matter fluid is inter- 
acting with the dark radiation component, the dark mat- 
ter component is no longer presureless, showing damped 
oscillations. The smaller wave mode at which the interac- 
tion between the dark fluids will leave a signature on the 
matter power spectrum is roughly kf afH{af), which 
corresponds to the size of the universe at the time that 
the dark radiation-dark matter interaction becomes inef- 
fective [H-il, i.e. when H{af) = r(a/) (being H the 
Hubble parameter and F the effective dark radiation-dark 
matter scattering rate -^j^ndm{o'd.m-dr\v\)) ■ For the case 
of constant dark radiation-dark matter interacting cross 
section, the typical scale kf reads, for ANgs = 1: 

^^/lO-32cm2MeV-i\'^\,, _i 

kf - 0.7 hMpc ^ , (13) 

\ Qq J 

Notice however from Fig. [T] (lower panel) that, while 
varying c^jg the matter power spectrum barely changes, 
a change in c^^ changes dramatically the matter power 
spectrum, washing out any interacting signature. For in- 
stance, if c^g = 0, dark radiation is a presureless fluid 
which behaves as dark matter, inducing an enhancement 
of the matter fluctuations, and, consequently, the pres- 
ence of a dark radiation-dark matter interaction will not 
modify the matter power spectrum, see the lower panel 
of Fig. [TJ Therefore, one might expect a degeneracy be- 
tween the dark radiation-dark matter coupling and the 
dark radiation c^^ parameter: the larger the interaction 
is, the smaller c^g should be to compensate the suppres- 
sion of power at scales k ^ kf. 

III. DATA 

In order to constrain the dark radiation parameters 
AA^cff, Cgjj and c^jg, as well as the size of the dark 
radiation-dark matter interaction, we have modified the 
Boltzmann CAMB code ^] including the dark radiation- 
dark matter interaction scenario. Then, we perform a 
Monte Carlo Markov Chain (MCMC) analysis based on 
the publicly available MCMC package cosmomc [s^. We 
consider a ACDM cosmology with AA'off dark radiation 
species interacting with the dark matter and three mass- 
less active neutrinos. This scenario is described by the 



4 



10000 



+ 



1000 



AN„„=1,c' 




=0.33, Ql 

=0, Qo=10"^^ err? MeV'^ 
WMAP7 datasets 
SPT datasets 



10 



100 



1000 



100000 



10000 



1000 



100 



10 



ANeff=1.c eff=0.33,c%i3=0.33,Q, 




LRG sample from Data Release 7 
ACDM WMAP7 

ANe„=1,cV0.Q( 



=0, Qo^ 



0.1 
0.0001 



0.001 



0.01 
k[h/Mpc] 



0.1 



FIG. 1: Upper panel: The magenta lines depict the CMB temperature power spectra C; for the best fit parameters for a 
ACDM model from the WMAP seven year data set. The dotted curve shows the scenario with a constant interacting cross 
section with Qo — 10~''^cm^ MeV~^ for AA^'cff = 1 and assuming canonical values for c^g = cJig = 1/3. The dashed (dot 
dashed) curve illustrates the same interacting scenario but with c^g — and c^is = 1/3 (c^g = 1/3 and Cvis = 0). We depict 
as well the data from the WMAP and SPT experiments, see text for details. Lower panel: matter power spectrum for the 
different models described in the upper panel. The data correspond to the clustering measurements of luminous red galaxies 
from SDSS II DR7 \M- 



following set of parameters: 

{Wb, Uc, &s, T, Us, log[10^°As], AiVcff, C^is, clg, Qo}, 

where ujb = ^th^ and Wc = Q,ch^ are the physical baryon 
and cold dark matter energy densities, 6s is the ratio be- 
tween the sound horizon and the angular diameter dis- 
tance at decoupling, r is the optical depth, Us is the 
scalar spectral index, is the amplitude of the primor- 
dial spectrum, ATVeff is the extra dark radiation compo- 
nent, c^ig is the viscosity parameter, c^g is the effective 
sound speed and Qq, in units of cm^ MeV~^, encodes the 
dark radiation-dark matter interaction. The flat priors 
considered on the different cosmological parameters are 



specified in Tab. ID 

For CMB data, we use the seven year WMAP data 
d, (temperature and polarization) by means of the 
likelihood supplied by the WMAP collaboration. We 
consider as well CMB temperature anisotropics from the 
SPT experiment [lo| . which provides highly accurate 
measurements on scales < 10 arcmin. We account as well 
for foreground contributions, adding the SZ amplitude 
Asz, the amplitude of the clustered point source contri- 
bution, Ac, and the amplitude of the Poisson distributed 
point source contribution, Ap, as nuisance parameters in 
the CMB data analysis. 

Furthermore, we include the latest constraint from 
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Parameter 


Prior 




0.005 ^ 0.1 




0.01 ^ 0.99 


e. 


0.5 ^ 10 


r 


0.01 ^ 0.8 


Us 


0.5 1.5 


In (10^" As) 


2.7 -5> 4 


AiVeft- 


^ 10 


Wis 


1 




1 


logio (Qo/lO-^'cm^MeV-i) 


-4^0 



TABLE I: Uniform priors for the cosmological parameters 
considered here. 



the Hubble Space Telescope (HST) [H on the Hub- 
ble parameter Hq. Separately, we also add Supernovae 
la luminosity distance data from the 3 year Supernova 
Legacy Survey (SNLS3) [H, adding in the MCMC anah 
ysis two extra nuisance parameters, which are related to 
the intrinsic supernova magnitude dependence on stretch 
(which measures the shape of the SN light curve) and 
color, see ReL jslj for details. We do not consider here 
the addition of HST and SNLS3 measurements simulta- 
neously because these two data sets are not independent. 

Galaxy clustering measurements are also added in our 
analysis via BAO data from the CMASS sample in Data 
Release 9 of the Baryon Oscillation Spectroscopic 
Survey (BOSS) [13, [ll], with a median redshift of z = 
0.57 [39], as well as from the LRG sample from Data 
Release 7 with a median redshift of z = 0.35 [1^, and 
from the 6dF Galaxy Survey 6dFGS at a lower redshift 
z = 0.106 tiU]. 

Therefore, we illustrate two cases, namely, the results 
from the combination of WMAP, SPT, SNLS3 and BAO 
data as well as the results arising from the combination 
of WMAP, SPT, HST and BAO data. 



IV. CURRENT CONSTRAINTS 

Table [H] shows the 68% and 95% CL errors on the 
dark radiation parameters and on the size of the dark 
radiation-dark matter interaction strength arising from 
the two possible combinations of data sets considered 
here for both interacting and non interacting scenarios. 
Notice, first, that the 1 — 2a preference found in the lit- 
erature for extra dark radiation species is still present in 
both interacting and non interacting scenarios in which 
the dark radiation clustering properties are not standard. 
Overall, the bounds on AA'cff are not largely modified 
when allowing for a dark radiation-dark matter coupling, 
see also the results presented in ReL [l8|. However, the 
bounds on the dark radiation clustering properties c^ff 
and in the ACDM scenario and in minimal extensions 
of this scheme presented in Refs. [H, |T3j are drastically 
changed when considering the possibility of an interac- 



tion between the dark radiation and the dark matter flu- 
ids. For instance, in ReL [l^, in the context of a ACDM 
scenario, it is found that c^g = 0.241°;°^ at 95% CL. Sim- 
ilar results were found in ReL flT^, where the ACDM sce- 
nario was extended to consider other cosmological models 
with a dark energy equation of state or with a running 
spectral index. Indeed, within non interacting scenarios, 
we find c^jf = 0.32;° °;^ and c^;^ = 0.271^-^2 at 95% CL 
from the combination of WMAP, SPT, HST and BAO 
data sets. These bounds are much weaker when allowing 
for an interacting dark radiation component: the errors 
on Cgff are degraded by an order of magnitude, while the 
errors on c^-^ increase by a factor of two. We find, for 
the same combination of data sets than the one quoted 
above, c^ff = 0.28t;^:28 and 0.45^^;^, both at 95% CL. 

Figure 111 (left panel) depicts the 68% and 95% CL ab 
lowed regions in the (c^fj, ATVcs) plane arising from the 
MCMC analysis of the cosmological data sets described 
in the previous section. We illustrate here the four cases 
shown in Tab. [TTl The green (yellow) contours refer to 
the case of WMAP, SPT, BAO and SNLS3 data sets 
with (without) interaction between the dark radiation 
and dark matter fluids. The magenta (red) contours re- 
fer to the case of WMAP, SPT, BAO and HST data 
sets with (without) interaction. Notice that the errors 
on the c^jf parameter are largely increased when the in- 
teraction term is switched on, while the errors on AN^s 
are mildly affected by the presence of such an interaction. 
Notice that HST data is more powerful than SNLS3 data 
in constraining AA^oft, agreeing with previous results in 
the literature, see ReL The reason is because AAes 
is highly degenerate with Hq, and HST data provide a 
strong prior on the former parameter. 

The right panel of Fig. H depicts the 68% and 95% CL 
allowed regions in the (c^jg, AN^g) plane, being the color 
code identical to the one used in the left panel. While 
the impact of the coupling is not as large as in the effec- 
tive velocity case, the errors on the viscosity parameter 
c^jg are enlarged by a factor of two in interacting dark 
radiation models. 



V. FORECASTS FROM FUTURE 
COSMOLOGICAL DATA 

We evaluate here the constraints on the dark radia- 
tion parameters, AA^oS, c^g, c^ig, by means of an anal- 
ysis of future mock CMB data for the ongoing Planck 
experiment and the future COrE mission. These CMB 
mock data sets are then fitted using a MCMC analysis 
to a non interacting cosmological scenario but allowing 
the dark radiation parameters to have non standard val- 
ues. The CMB mock data sets are generated accord- 
ingly to noise properties consistent with the Planck and 
COrE CMB missions. The fiducial Cg model we use is 
a ACDM scenario (i.e. a flat universe with a cosmologi- 
cal constant and three massless active neutrino species) 
adding an interaction between the dark radiation and 
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WMAP+SPT+BAO2012 WMAP+SPT+BAO2012 WMAP+SPT+BAO2012 WMAP+SPT+BAO2012 
+HST int. +HST no int. +SNLS3 int. +SNLS3 no int. 



^2 

Ccff 



0.28 



+0.10+0.44 



0.32 



+0.02+0.04 



0.30 



+0.12+0.50 



0.32 



+0.02+0.04 



0.27 



+0.13+0.34 



68%CL 
95%CL 



< 0.81 

< 1.30 



0.62 



+0.36+0.80 



< 0.76 

< 1.47 



r> 77+0.50+1.29 
' ' -0.72-0.72 



Qo 



68%CL 



(IQ-^^cmVMeV"^) 95%CL 



< 0.8 

< 4.9 



< 0.8 

< 5.4 



TABLE II: ID marginalized bounds on the dark radiation parameters and on the size of the dark radiation dark matter 
interaction Qo using WMAP, SPT, BAO data and IIST/SNLS3 measurements, see text for details. We show the constraints 
for both interacting and non interacting models, presenting the mean as well as the 68% and 95% CL errors of the posterior 
distribution. 




FIG. 2: Left panel: 68% and 95% CL contours in the (c^g , AiVeff) plane arising from the MCMC analysis of WMAP, SPT, 
BAO and HST/SNLS3 data. The green (yellow) contours refer to the case of WMAP, SPT, BAO and SNLS3 data sets with 
(without) interaction between the dark radiation and dark matter fluids. The magenta (red) contours refer to the case of 
WMAP, SPT, BAO and HST data sets with (without) interaction between the dark radiation and dark matter sectors. Right 
panel: as in the left panel but in the (Cvis, AA'^off) plane. 



dark matter sectors with Qo = 10~^^cni^ MeV""'^, as- 
suming one dark radiation interacting species AiVeff=l 
and standard clustering and viscosity parameters for the 
dark radiation, i.e. c^jg = clg = 1/3. For each fre- 
quency channel we consider a detector noise given by 
a;-i = {eaf, being 6 the FWHM of the gaussian beam 
and a = AT/T the temperature sensitivity (the polariza- 
tion sensitivity is AE/E = y^AT/T). Consequently the 
Ce fiducial spectra get a noise contribution which reads 
Ne = cj-i exp (£(£ + where h = 

Figure [3] (left panel) depicts the 68% and 95% CL con- 
tours in the {clg,ANes) plane arising from the MCMC 
analysis of Planck and COrE mock data. Notice that the 
reconstructed value for c^g is larger than the simulated 
value clff = 1/3. The reason for that is due to the degen- 
eracy between the dark radiation-dark matter interaction 



Qo and c^g, see Fig. [1] from which one would expect a 
negative correlation between the interaction cross section 
and the effective velocity. If such an interaction occurs in 
nature but future CMB data is analysed assuming a non 
interacting model, the reconstructed value of c^g will be 
higher than the standard expectation of 1/3, see Tab. lIIII 
From what regards to c^j^, see Fig. [3] (right panel), the 
effect is the opposite since these two parameters are posi- 
tively correlated and therefore the reconstructed value of 
c^ig is lower than the canonical 1/3, see Tab. IIIII There- 
fore, if the dark radiation and dark matter sectors inter- 
act, a large bias on the dark radiation clustering param- 
eters could be induced if future CMB data are analysed 
neglecting such coupling. On the other hand, the bias 
induced in AA^cff is not very significant, being the recon- 
structed value consistent with the AA^eff = 1 simulated 
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one within la. 



VI. CONCLUSIONS 

Standard dark radiation is made of three light active 
neutrinos. However, many extensions of the standard 
model of elementary particles predict an extra dark radi- 
ation component in the form of sterile neutrinos, axions 
or other very light degrees of freedom which may interact 
with the dark matter sector. In fact, once that one as- 
sumes the existence of extra dark radiation species as well 
as the existence of a dark matter sector there is a priori no 
fundamental symmetry which forbids couplings between 
these two dark fluids. If one allows for such a possibil- 
ity, the clustering properties of these extra dark radiation 
particles might not be identical to those of the standard 
model neutrinos, since the extra dark radiation particles 
are coupled to the dark matter. In this paper we have an- 
alyzed the constraints from recent cosmological data on 
the dark radiation abundances, effective velocities and 
viscosity parameters. While the bounds on ANcS are 
very close to those of uncoupled models, the errors on the 
clustering dark radiation properties are largely increased, 
mostly due to the existing degeneracies among the dark 



radiation-dark matter coupling and c^fj, cjjg. The cos- 
mological bounds on the dark radiation effective veloc- 
ity Cgfj found in non-interacting schemes are degraded 
by an order of magnitude when a dark radiation-dark 
matter interaction is switched on. In the case of the vis- 
cosity parameter c^jg, the errors on this parameter are 
a factor of two larger when considering interacting sce- 
narios. We have also explored the perspectives from fu- 
ture Cosmic Microwave Background data. If dark radi- 
ation and dark matter interact in nature, but the data 
are analysed assuming the standard, non interacting pic- 
ture, the reconstructed values for the effective velocity 
and for the viscosity parameter will be shifted from their 
standard 1/3 expectation, namely c^g = 0.34to;oo3 and 
c2.^ = 0.29t°;°°2 at 95% CL for the future COrE CMB 
mission. 
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FIG. 3: Left panel: 68% and 95% CL contours in the (Cgff , AA'eff) arising from the MCMC analysis of Planck (red contours) 
and COrE (blue contours) CMB mock data. The mock data are generated adding an interaction between the dark radiation 
and dark matter sectors with Qo ~ 10~^^cm^ MeV~^, assuming one dark radiation interacting species AA'cff^l and standard 
clustering and viscosity parameters for the dark radiation. The CMB mock data is then fitted to a non interacting cosmology 
but allowing the dark radiation parameters c^g and cJi^ to have non standard values. Right panel: as in the left panel but in 
the (c^is, AA'eff) plane. 
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